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The development of the brush seal is considered to be most promising amongst the 
advanced type seals that axe presently in use in the high speed turbomachineiy. The 
brush is usually mounted on the stationary portions of the engine and has direct 
contact with the rotating element, in the process of limiting the "unwanted” leakage 
flows between stages, or various engine cavities. This type of sealing technology is 
providing high(in comparison with conventionaLseals) pressure drops due mainly to the 
g packing density(around 100 bnstles/1 mm 2 ), and brush compliance with the rotor 
motions. In the design of modern aerospace turbomachinery leakage flows between the 
stages must be minimal thus contributing to the higher efficiency of the engine. Use of 
the brush seal instead of the labyrinth seal reduces the leakage flow by one order of 
magnitude^, 2], Brush seals also have been found to enhance dynamic performance, cost 
less and are lighter than labyrinth seals. Even though industrial brush seals have been 
successfully developed through extensive experimentationfl,2l, there is no 
comprehensive numerical methodology for the design or prediction of their 
' Tte ®? astm S analytical/numerical approaches are based on bulk flow 
models[6,7J and do not allow the investigation of the effects of brush morphology f bristle 
arrangement) or brushes arrangement (number of brushes, spacing between them) on 
the pressure drops and flow leakage. An increase in the brush seal efficiency is clearly a 
complex problem that is closely related to the brush geometry and arrangement, and 
can be solved most likely only by means of a numerically distributed model. 
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STATE-OF-THE-ART 


★ The reduced leakage, and physical compliance of the brush 

body to external perturbing factors are features that stand 
out in turbomachinery applications where there are f 

expected boundary variations due to mass flow, brush fibers 
compliance pressure, temperature, and time dependent 
eccentric shaft motion. All these characteristics have made the 
brush configuration an especially interesting and worthy 
candidate. 

★ Rolls-Royce(RR), in 1980’s, has successfully introduced a brush 
seal manufactured by Cross Mfg. Ltd.(CML) on a 
demonstrator engine, and then tested it for several thousand 
hours, Fergussonfl]. More recently EG&G Sealol, Technetics, 
Detroit-Allison and others have enabled full programs of study 
of this type of seal. 

★ Conclusions of a recent workshop on code development(1992) 
indicate that while the brush seals works well, there is a need to 
improve its performance characteristics. Such a goal can be 
achieved by using cascades of brushes, nonhomogeneous brush 
morphology, "non-conventional" brush structure design, and in 
general, a process of optimization of brush design parameters. 

★ The concept employed by the lumped bulk flow numerical 
models can not predict local brush compliance, associated local 
flow phenomena and the pressure drops and the transient 
effects associated with them. The importance of the local flow 
phenomena in the sealing process is paramount to the global 
performance of the brush. 
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CURRENT RESEARCH ACTIVITIES : 


• DEVELOPMENT AND VALIDATION OF A NUMERICAL 
ALGORITHM AND COMPUTER CODE THAT UTILIZES 
MATHEMATICAL MODEL WITH DISTRIBUTED PRIMARY 
PARAMETERS(NAVIER-STOKES EQUATIONS), 

UNDER NASA GRANT , NASA LEWIS RESEARCH 
CENTER. 


Computer Code Allows: Estimation of the pressure drops (flow 
rates) for the typical brush seal segments of different shapes: 
i.e. brustles diameters, configurations, packing densities. 


APPROACHES: 

-Large size characteristic segment: 7-10 rows with 10 pins in one 
row in the transversal direction 

-Brush Partitioning: inflow segment(first 3 rows), central part, 

outflow segment(last rows) 


• DEVELOPMENT OF THE COMPUTER PROGRAM THAT 
ADRESSES BRISTLES MOTION AND ITS INFLUENCE ON 
THE PRESSURE DROPS AND FLOW RATES. 


• FURTHER MODIFICATION OF THE EXPERIMENTAL 
FACILITY FOR THE PURPOSES OF CODE VALIDATION 
DIFFERENT SHAPES OF THE BRUSH SECTION 
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OBJECTIVES AND ACCOMPLISHMENTS: 


.Develop verified family of CFD codes for Analyzing Brush Seals 

-idealized(uncompliant) 2D configuration \/ 

regular gridding \/ 

variable grid size \/ 

-compliant 2D geometry under development 

• Experimental Facilities for the Adequate Code Verification 

-stationary bristles(cylinders) \/ 

-moving bristles under development 

• Qualitative and quantitative analyses of the Fluid Flow in 
the Brush Seal Configuration 

-flow around one bristle, level 1 \/ 

-flow around several bristles, level 2 \/ 

-flow in deep tube bundles, level 3 \/ 

(intermediate pitch-to-diameter ratio) 

7 rows of pins with 11 cyl. in a row 


-flow through uncompliant brush prototypes 
(small pitch-to-diameter ratio), level 4 \/ 

-flow through the characteristic brush 
segments, brush partitioning technique \/ 

level 5 
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PROBLEM DEFINITION 


Conclusions of a recent NASA Seal Workshop[5] indicate that while 
the brush seals work well, there is a need to further improve the 
performance characteristics. Such a goal can be achieved by using 
cascades of brushes, nonhomogeneous brush morphology, "non- 
conventional” brush structure design, and in general, a process of 
optimization of brush design parameters[18]. The distributed 
velocity fields(u,v) and the associated pressure maps are of vital 
importance for the prediction of the average pressure drop, or the 
possible sudden failure of the brush seal under unexpected local 
pressure hikes”. The momentum carried by these velocities(or the 
upstream pressure) can force the brush deformation, and can 
create favorable conditions for the brush 'opening', followed by seal 
failure[4,15]. It is in this context that the development and 
validation of a numerical model with distributed primary 
parameters(u,v,p) becomes important. 


The design goals of the model are to determine the pressure drop 
for a configuration specified by the designer, i.e. the density of the 
brush packaging, length, number of rows, bristles 
sizes(homogeneous or nonhomogeneous), distances between the 
rows, or brushes respectively(single, double or cascade brush). The 
systemic goals are a) to develop physically relevant packages of 
assumptions for the simulation of the brush seal and to implement 
a robust numerical methodising primitive variables u,v,p ) for the 
calculation of the forced convective flow through dense brush-like 
cylinder arrays, and b) to analyze numerically different aspects of 
the flow dynamics in the generic brush prototypes. Achieving the 
goals set forth in a) and b), will allow usage of predictive design 
codes with a high level of reliability. 

For analysis and classification purposes we identified four 
different models of the generic brush geometry[19]. The simplest 
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model (Level 1) assumes flow analysis in the vicinity of a non- 
moving single brush bristle[20]. The second level(Level 2) model 
introduces the analysis of a limited cluster that consists of several 
non-moving bristles[19,20], that do influence jointly the flow field 
through disturbances generated by the wake vortices. The next 
level(Level 3) introduces the analysis for a multi-cluster[21] f and 
finally we assemble large numbers of clusters with small pitch to 
diameter ratio(PTDR) that actually simulate the real brush( Level 4, 
[22] ). Each one of these levels is designed to introduce one 
additional level of difficulty, help learn more about the physics of 
the flow, and increase the level of confidence in the final numerical 
model. 

One can see an idealized schematic of a linear brush seal in 
Fig. 1. In real conditions the flow upstream of the seal exhibits 
both circumferential and axial components. Reynolds numbers that 
are typical for the circumferential component are usually in the 
range 10 4 - 10 6 , while the Reynolds number of the transversal 
component component (leakage flow) does not exceed low laminar 
values that are defined by the design requirements of the seal (an 
ideal case leakage is equal to zero). A review of the data published 
by Chupp et al[6], Nelson and Chupp [10], Dowler[ll] allowed us 
to determine the typical ranges of the parameters that are usual for 
the brush seals tested by the industry. The level of the leakage flow 
and maximum velocity in the pitch between the bristles were 
estimated as 3.17 m/s and 76.55m/s respectively. In Fig. 2 one 
can find an approximate range of Reynolds numbers typical for a 
brush seal functioning in air. As we can see from this qualitative 
figure an assumption of laminar and incompressible fluid can be well 
justified since the Reynolds numbers are not in the turbulent range 
and the Mach number is M<0.4. The authors have established 
through their experimental work[4,15] that the major factor 
contributing to the pressure drop is the longitudinal flow in the X 
direction in the XOY plane(Fig. 1). Thus, in order to simplify the 
model we assumed a two dimensional flow, and neglected both 
curvature effects and flow in the Z direction. 
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Figure l(a,b) Idealized Schematic of the Brush Seal and Flowpath 

a) generic seal 

b) problem definition 
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M = 0.3 M = 0.4 



Figure 2. 


Qualitative Diagram of the Flow Leakage vs 


Reynolds Number 
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MATHEMATICAL MODEL 


• Characteristic Geometry 

The computational domain (with L x as length and L y as width) is 
represented by the horizontal brush cross section restricted by two 
walls as shown in Fig. lb. Inside the domain, solid bodies(bristles) 
of round cross-section are located, thus, creating a structure of (n) 
rows with (m) elements in each row. The bristle diameter is used as 
a characteristic length scale(LQ) and the velocity at the entrance, 
as a characteristic velocity(UQ). 

• Governing equations. 


It is assumed that the flow is viscous and laminar and it is caused 
by the pressure difference across the generic brush. The initial 
velocity distribution u(x,y) and the characteristic pressure at the 
entrance are assumed known. The influence of the body forces is 
excluded from the problem definition. We also neglect the 
influence of heat transfer on the flow structure(isothermal flow). 
The two-dimensional Navier-Stokes equations for unsteady 
incompressible viscous flow, can be written in dimensionless 
conservative form (Cartesian), as 


au , a(m*) , d(uy) _ 
at ^ ax + ay ~ 


a P 
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(3) 

(4) 



MATHEMATICAL MODEL 


During calculations, -s(the dilation term) represents the residual 
that has to shrink to zero if continuity is to be satisfied[23,24]. 
Previous numerical work of these authors[19,20,21] exp ore 
models of Levels 1, 2, and 3 in order to obtain qualitative 
information about fluid flow in the generic brush configurations. 

The new technique that involves brush partitioning proposes the 
division of the brush in three typical areas(Fig.l): (i)inflow area 
with the free flow and first several rows of pins with developing 
flow profile, (ii) a central core that consists of several rows of 
bristles with symmetric fully-developed flow distribution, and (iii) an 
outflow zone that includes the last rows of bristles and the area 
behind the brush where vortices generation is taking place. 

Hendricks[25l showed that the first thre e and^ last thre e ro^ 
bristles are usually deflected outward of the main brush core during 
the experiments. These bristles are located at a considerable 
distance from each other and thus do not significantly contribute to 
the overall pressure drop. Since the central core has the highest 
level of packing density and flow resistance, as a first step we wi 
concentrate our analysis on this region. In general, the tota 
pressure drop in cylinder arrays can be expressed as following: 

aP— C z ^APgy , i— — 1 - * ■■■■* 

According to Zukauskas et al.[18], C =1. if the number of rows is 
greater than three. One can conclude that in terms of the pressure 
field the flow is fully developed as it passed the first three rows. 

Assumption of a fully developed flow distribution in the central 
brush core reduces the size of the computational domain to a 
characteristic cell of several sequential rows(characterized by a 
constant aP 3 v). if indeed, one supplies proper symmetric 
conservative boundary conditions. 
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BOUNDARY CONDITIONS 


• Brush Entrance Region 

We assume an entrance velocity profile with u=/(y,t) and v=/(y,t) f 
which in most cases can be reduced to u(x=0)=l and v(x=0)=0. 
In the case of lateral solid walls boundary conditions are specified 
as non-porous and non-slip, i.e. u=v=0. In the case of flow 
symmetry, boundary conditions at the "top" and ” bottom" 
boundaries are given as 

a) fjf=0 and ^=0 (6) 

b) ^=0 and v=0 (7) 

Exit boundary conditions can be specified in a similar way: 

§X = ® a °d v=0 (8) 


• Brush Core Zone 

Inflow velocity condition for this zone are the outflow conditions for 
the brush entrance zone. In this case: 

v=0 and u^fy.t), (9) 

where fi(y,t) is established from the solution of the Navier-Stokes 
equations in the entrance region. At the outflow of the 
computational domain we apply Eq. 8, and at the lateral walls Eqs. 
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boundary conditions 


• Brush Outflow Zone 

Inflow velocity conditions for this area are calculated via the 
solution of the Navier-Stokes equations at the brush core zone. 
That implies v=0 and u=f 2 (y,t). At the lateral walls we employ 
conditions (6-7). The exit flow is allowed to develop naturally based 
on boundary conditions situated far downstream The velocity 
boundary conditions imposed at the exit are derived from (i) t e 
satisfaction of the continuity equation(specifically in the direction ot 

the U velocity), 

au _ ay (10a) 

ax — ay v 7 

and (ii) the condition for fully developed velocity gradient for the 
vertical component 


s =0 ( 10b ) 

• Boundary Conditions for the Poisson Equation 


The dynamic pressure P, is determined from the balance of the 
normal forces with the inertia and viscous forces. This formulation 
implies Neumann type boundary conditions. The effects of e 
terms ^ and ^ have been considered negligible. The following 
formulation o^Egs. 11 and 12 is totally independent of the 
boundary configuration, and thus applicable with no restrictions. 


aP = _ I 

ax 1 


aP _ f 

ay — 1 


a(uu) 

ax 


a(uv) 

ax 


+ 
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a(uw) 1 , 9 2 u , a 2 u y 
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a(w) 1 ( 0 2 v + » 2 X) 
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In combination with the Dirichlet conditions for the velocity, Eqs. 
11 and 12 form a well-posed boundary-value problem with resulting 
solutions free of odd-even velocity-pressure decoupling as described 
by Gresho[26] and Patankar[27]. 


Odd-even Velocity- Pressure Decoupling For the Inline Brush 

Segment 



Initial conditions(t=Q). The input velocities are given as u=l, v=0. 
The pressures are set initially to an arbitrary, operator chosen 
P=P re f=const. 
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SOLUTION PROCEDURE 


The discretization of the system of governing equations introduced 
above, follows the use of the Alternating Direction Method[28| 
applied to a collocated grid. The procedure uses the full direct 
approximation of each term within the differential equation on 
every half time step, a,/2. One obtains the following system of 
linear algebraic equations. 

The spatial derivatives, with the exception of the convection terms 
and cross-derivatives, are approximated by an implicit second order 
central finite difference. For the convective terms, the implicit form 
of the third-order deferred correction scheme proposed by 
Kudriavtsev[30] and Hayase et al.[3l] was implemented[4]. 


After discretization of the boundary conditions we obtain a self 
sufficient system of linear algebraic equations that is solved using a 
tridiagonal matrix elimination. The steps of the solution are as 


• introduce the initial flow and pressure field at the n time level 

• solve in the x direction for U velocity at the n+1/2 time step 
. solve in the y direction for the U velocity at the n+1 time step 

. solve in the x direction for V velocity at the n+1/2 time step 
.solve in the y direction for the V velocity at the n+1 time step 
. solve the Poisson's pressure equation at the n+1 time step by 
means of the pseudo-transient method within the set of internal 
ADI iterations n it =l,...n f : (i) in the x direction at the s+1/2 

pseudo-time step and (ii) in the y direction at the s+1 pseudo- 


time step 

• advance to the next time level . 


On each iteration pressure at the reference point P(l,2) is assigned 
as P*. Nondimensional pressure at every node is calculated as. 
A p . _ p* _ p_ |f steady-state solution is of interest we employ 

psevdo-transient method as proposed by [28]. 
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RESULTS AND DISCUSSION 


• Pressure and flow patterns for different generic brush formations. 

There is a considerable amount of numerical and experimental work 
treating the flow around circular or squared cylinders in crossflow. 
However, only few investigations concentrated their attention on 
the flow interaction with the tube bundles[32,33,34J. The effects of 
the cylinders' arrangement and the array size, or morphology on 
the flow structure and pressure drops has not been extensively 

studied, especially if the pitch-to-diameter ratio(PTDR) is smaller 
than 1. 

• Braun et al.[21], Kudriavtsev et al.[20], and Kudriavtsev[19] have 
studied systematically the time-flow, and time-pressure 
development in arrays of cylinders with PTDR equal and smaller 
than 1. The database that emerged, showed that the development 
of . flow and pressures around one cylinder, small groups of 
cylinders, and large groups of cylinders, do not lend themselves to 
extrapolation from one configuration to another. This is exemplified 
in Figs. 3 and 4. Figure 3 presents the pressure oscillations at the 
trailing edge of a single cylinder located in the square channel of 
Fig. 1. On inspection, it can be seen that in a quasi steady regime 
the. pressures in the wake of the cylinder reach a repetitive 
oscillatory profile within an envelope of +5%. When the number of 
cylinders is increased to five, arranged in a staggered formation of 
two rows, the same envelope is contained between +10%, Fig. 4a. 
Finally, the increase of the cluster to 72 cylinders, arranged also in 
a staggered formation of seven rows, generates in the wake of the 
array oscillatory pressures that can be contained within +2.5% 
envelope, fig. 4b. These results, when considered in conjunction 
with the flow patterns already discussed by the authors in previous 
papers, demonstrate that phenomena indigenous to small arrays do 
not, as a rule, apply to large formations. Therefore one has to be 
extremely careful during the analytical/numerical modeling of large 
arrays. Special methodologies, that involve realistic treatment of 
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PARAMETRIC STUDIES OF THE FLOW IN THE 
BRUSH PROTOTYPES LOCATED IN THE OIL TUNNEL 

Re= 






I 




/a.1 


■got 




grid 551x1281 qrid 501x1281 

CALCULATED FLOWFIELD 
tP(num) =4.31 psi AP(exp)=4.02 psi , error =7% 









non dimensional time, t 
b) 

Figure 3(a,b) Pressure Development Behind Single Cylinder(Level 1) 
and Small Cluster of Cylinders (Level 2) 

a) single sylinder in the square channel 

b) small cluster in the square channel 
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boundary conditions, have to be applied to small arrays, that can 
then be assembled (based on these boundary conditions) in large 
clusters, and render a realistic reproduction of the flow conditions. 

Figures 5 present the flow development in an array of 
7(rows)xll round pins with PTDR^=PTDRy=l. Along the walls, 
one can clearly observe high velocity rivering jets that engender 
contiguous regions of lower pressure which attract the flow from 
the center of the generic brush[21]. The central region appears to 
be, fortunately, quite repetitive in nature, and indicates that the 
use of a repetitive cell for the construction of a large brush core is 
a feasible alternative. Finally the wake region presents the 
formation of the exit jets, and the incipient formation of 
meandering vortices. 

The same algorithm used to obtain the results presented in 
Fig. 5(PTDR=1) has been used for the modeling of the dense 
array of Fig. 6. The flow, at Re=195, presents a steady-state 
pattern in this array of three rows with PTDR^=PTDRj=0.084. 
This packing density is of significance for the real industrial brush 
seal applications(l/10>PTDR>l/30). Modeling of the fluid flow 
through the pin array with such a PTDR is a challenging 
computational task that so far has not been addressed in the open 
computational fluid dynamics(CFD) literature. The flow is 
dominated by accelerating-decelerating streams formed by the 
system of converging-diverging nozzles(constituted by the adjacent 
cylinders). These flows are characterized by extremely high 
convective gradients at their smallest cross-section. This situation 
imposes limitations on the computational algorithm in terms of 
computational stability, resolution and global accuracy. The window 
marked in Fig. 6a, is shown in detail in Fig. 6b. One can get an 
appreciation of the complicated flow formations that are at work. 
In Fig. 6b the rivering, lateral rivering, and the nozzle-jets are quite 
evident. The formation of the jets at the exit of the array can also 
be observed. The corresponding experimental results and the 
successful superposition of these results and the numerical 
simulation are shown in this sequence(Fig. 6b). For an entire array 



node number 



Computational Grid 139x139 


Figure 5. Typical Flow Structure in the Large Cluster of Cylinders(Level 3) 
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of 3(rows)xll, Fig. 7a, presents experimentally the interaction of 
12 convective jets as they interact with each other and with the 
wall. The experiment was run at Re=1000, and a rather unusual 
flow formation was observed in the wake of the array. The low 
pressure region that developed near the walls, in- and outside of 
the pin array, causes constant jet flow deviation towards the walls. 
The void created in the wake is associated with a pressure that is 
lower than the far downstream pressure. This pressure field 
engenders a recirculation bubble that occupies considerable part of 
the outflow cross section and consists of the two butterfly like 
vortices with strong "negative” flow along the centerline of the 
cross section. The global resulting recirculation structures include 
numerous vortices and require high computational accuracy in order 
to obtain any realistic solution. Such solution can be observed in 
Figs. 7b and 7c. For the numerical results of both Figs. 6 and 7 it 
has been found that the flow structure, inside the array, converged 
to a steady solution much faster than the flow in the wake of the 
array, Figs. 7b, 7c. 

For PTDR<0.1, the numerical solution can be obtained only 
by large scale CFD modeling that requires grid sizes approaching 
1000x1000 nodes. Thus a further increase in the size of the array 
or any further grid refinement becomes rather prohibitive, and 
appropriate strategies need to be found to handle this situation. 
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b)numerical(smaU array) 



c)numerical(large array) 

Figure 7(a,b,c) Butterfly Formations in the Wake of the Cylinder’s Clusters 

Re=1000 
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PARAMETRIC STUDIES OF THE FLOW IN AN ARRAY OF 
FIVE STAGGERED CYLINDERS 


authors have extended their analysis to a representative cell of 
ive staggered cylinders. The longitudinal and transversal pitches 
are equal to each other, PTDR l =PTDR t = 1, while the cylinder 
diameter is d=0.143L n . 



CORRESPONDING 
BUTTERFLY 
FLOW FORMATION 


A 



LAST ROW OF PINS 
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The new approach of brush segmentation and assembly 


The totality of the numerical and experimental work presented 
above has led us to the conclusion[19], that we were fastly 
approaching the limits of our modeling capability from the 
standpoint of the computer capacity to handle the large arrays 
necessary to model a full brush. We have learned that 

a) the flow inside the brush reaches a convergent solution, much 
sooner than near the wall boundaries or in the brush wake; 

b) in a deep cluster the flow is practically repetitive, as long as it is 
at least two rows away from the inlet or exit of the generic brush; 

c) the pressures fall monotonously in the core portion of the 
simulated brush[4]; 

d) the computational time for calculating flow in dense arrays of 
cylinders(3rows x 11 lines) at moderate Re numbers can become 

prohibitive. 

As mentioned earlier, the brush can be partitioned into several 
areas. The numerical simulations shown in Figs. 8 explain the 
rationale of the brush partitioning approach. One of our previous 
conclusions was that large brush modeling is restrictive, so one tries 
to reduce the computational domain to a representative cell. In 
order to define optimally this cell, and its environment two 
Geometrical setups were chosen. First is the flow in the channel 
with non-slip walls as boundary conditions. The second setup 
replaces the wall with symmetry conditions at the boundary. From 
Fig.8(a,b) one can see that the wall introduces near wall effects 
that considerably change the flow structure. Thus one has to select 
the symmetry boundary conditions setup, as the optimal one(see 
Eqs 7-8 10) At Re=100 in the later time stages of the flow 
development one observes the incipience of an asymmetric mode in 
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the outflow zone. This situation induces asymmetry in the central 
part of this brush segment(Fig. 8b and 8d). In larger(deeper) 
cylinders bundles this effect does not penetrate to the core(like in 
Fig. 5). Thus, the core and the outflow segments need to be 
treated as separate computational domains. 

To respond to these physical conditions and the conclusions 
drawn above in items a) to d) we proceeded to assemble a full 
brush out of the component segments that were identified in the 
mathematical model section. The continuity and momentum 
governing equations remain the same, but the boundary conditions 
had to be modified in order to take into the account the repetitive 
nature of the assembly of the brush core. These boundary 
conditions were described by Eqs. 6-10. For the verification of the 
concept we have used a partitioned array that contains the 
segments shown in Fig. 9. Figure 9a, presents a succession of 
identical core segments that are characterized by the same pressure 
drops. Figure 9b displays the inflow and outflow segments that are 
added to the central core in order to form a complete system. The 
resulting total pressure drop can be calculated as 

aP brush = aP inflow + N** 9 aP central + aP outflow ( 19a ) 

If one considers the Hendricks[25] report, and if the brush has 
more than n>io rows, it seems rational that the overall pressure 
drop can be evaluated with little error by considering only what we 
defined as the core portion. Thus one can simplify Eq. 19a to 

AP brush = n aP centra | (19b) 

where AP * cen t ra | ~ average pressure drop per single row. 
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Figure 8(a,b,c,d) Samples of the Fluid Flow in the Brush Segments at Re-100 

a) symmetric mode(uutial stage of flow development) 

b) asymmetric mode(final stage of flow development) 

c) symmetric mode(initial stage of flow development) 

d) asymmetric mode(final stage of flow development) 







STUDIES OF THE FLOW IN THE BRUSH SEGMENTS 


Variation in the boundary conditions 



Re d =20 

Boundaries are opened to the flow 
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STUDIES OF THE FLOW IN THE BRUSH SEGMENTS 


boundary pitch equal to the internal pitch 




assymetric configuration, top pitch larger than the bottom 

Re_i=20 



Cooling Pins : Stages of the Transient Flow Development 
OTHER POTENTIAL APPLICATIONS OF THE CODE: 



, 1^=3000 (dimensionless time T=h^/UQ*t —0.21 ) 
fte H =1000 

U o =0.5 m/s ,» 0 =995 (working fluid - oil) 
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CONCLUSIONS 


The authors have presented chronologically the main contributors 
in the technological development of the brush seals. A 
see there were three avenues that were followed The first involved 
laboratory experiments on simulated brush seals, that led to 
better understanding of the flow conditions and flow paths in the 
brush The second involved industrial testing of brus , 
yielded real experimental data concerning pressure drops i and 
Lakaaes These efforts performed within the framework °f major 
industrial manufacturers and users, led to the incorporation of these 
S in production jet engines. The third and final avenue was the 
numerical development of design and predictive t 
brushes These codes have used either lumped or distributed c , 
« ft was shown earlier in this paper. The present paper presents a 
first attempt at extending the distributed model of Braun and 
Kudriavtsev to a model where the brush is segmented in inlet, core 
and exit component. This concept allows the construction of 
infinitely larqe brushes with diminished computational penalty It 

IS*. *>■ IM. concept p.«fc.s cnccll, ^JT.S 5 
repetitive segments of the core, as well as m the inlet and exit 

zones. 
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OBTAINED NEW SCIENTIFIC RESULTS : 


* transients are negligible for Re < 1000 

* pressure distribution and flow structure behind the simulated 
brush(located in the channel) are non-symmetrical. Flow 
assymetry depends on the Re number. 

* creation of the butterfly-like flow formation behind the simulated 
brush seal(located in the channel) and observed experimentally 

and numerically 

* nonlinear behaviour of the pressure coefficient vs Re number 

for the brush segments located in the channel. Pressure paradox. 

* nonlinear behaviour of the pressure drop vs PTDR L 

* formation of the near wall jets for the multi-row cylinder bundles 
located within the channel (simulated brush), Coanda effect, flow 
"expansion” from the center to the walls 


196 


CONCLUSIONS 


★ approximate mathematical model with distributed parameters 
that is free from the disadvantages associated with porous media 
assumption 


★ developed and evaluated a new computational algorithm for the 
solution of the N-S equations in (u,v,p) formulation 


★ systematic analyses of the fluid flow in the brush seal 
components 


★ Capabilities of calculation of the pressure drop for a given bristle 
geometry and the flow rate for typical brush segment 
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